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Nulling  in  Linear  Array  Patterns  With 
Minimization  of  Weight  Perturbations 


1.  INTRODUCTION 

Adaptive  antenna  systems  presently  receive  much  Interest  and  there  has  been 
considerable  study  of  adaptive  algorithms.  Less  attention  has  been  devoted  to  the 
end  product  of  the  algorithms,  which  is  an  antenna  pattern  with  nulls  In  the  direc¬ 
tions  of  the  jammers.  In  this  report  we  therefore  study  the  general  properties  of 
patterns  with  null  constraints.  The  hope  is  that  this  may  contribute  to  more  insight 
into  the  basic  behavior  and  fundamental  limitations  of  adaptive  nulling  systems  and 
possibly  also  suggest  new  adaptive  schemes  suitable  for  large  array  antennas. 

The  approach  we  take  is  based  on  a  direct  synthesis  of  the  array  pattern  with 
the  desired  nulls,  in  contrast  to  an  iterative  solution  using  an  adaptive  algorithm. 
The  problem  is  formulated  as  an  approximation  problem:  a  set  of  complex  array 
weights  JwnJ  Is  determined  which,  subject  to  the  null  constraints,  best  approxi¬ 
mates  a  given  set  of  weights  |won}  •  Two  types  of  approximation  criteria  are 
considered:  (a)  minimum  relative  weight  perturbations;  that  is, 

T.  |(w_  -  w  )/w  | 2  »  mir. ,  and  (b)  minimum  total  weight  perturbations;  that  is, 

n  on  on 

£1  wn  -  wQnl 2  *  min.  Criteria  related  to  (b)  have  been  applied  before, 2  whereas 
(Received  for  publication  18  February  1982) ' 

1.  Mayhan,  J.  (1976)  Nulling  limitations  for  a  multiple-beam  antenna,  IEEE  Trans. 

Antennas  Proprg.  AP  -24:769-779. 

2.  Steyskal,  H.  (1982)Synthesis  of  antenna  patterns  with  nulls,  IEEE  Trans. 

Ar.’ennas  Propas.  AP-30, 

. — ,i  —  ■  ■—■*■■■»  W 
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the  criterion  (a)  teems  to  be  new.  For  both  cates  exact  solutions  are  derived  and 
interpreted  in  terms  of  cancellation  beams  superimposed  on  the  initial,  uncon* 
strained  pattern. 

We  then  consider  the  same  synthesis  problem  under  the  restriction  that  the 
array  excitation  is  perturbed  in  phase  only.  The  motivation  for  this  is,  of  course, 
that  in  a  phased  array  the  required  electronic  controls  are  available  at  no  extra 
cost.  This  subject  is  of  great  practical  Interest,  and  several  publications  on  it 
have  appeared  in  the  literature.  The  associated  approximation  problem  is  non¬ 
linear  in  general  but  can  be  linearised  by  assuming  that  the  required  phase  perturba¬ 
tions  are  small,  an  assumption  which  is  reasonable  in  certain  applications  involving 
null  placement  in  low  sldelobe  regions.  Based  on  this  assumption  we  determine  a 
set  of  phase  perturbations  j  pn|  for  the  same  two  criteria  as  above  which  now  take 

the  forms  (a)  minimum  phase  perturbations:  that  is,  =  min. ,  and  (b)  minimum 

products  of  the  phase  perturbations  with  the  element  Amplitudes;  that  is, 

»I-J^  n  min.  This  time  the  solutions  cm  be  interpreted  as  pairs  of  can¬ 
cellation  beams  superimposed  on  the  initial  pattern. 

Following  the  analysis  of  these  null  synthesis  problems  we  present  the  results 
of  comi  /<ations  designed  to  display  the  basic  features  of  the  solutions  obtained  and 
to  explore  the  limitations  of  the  small  phase  perturbation  assumption  used  to  linear¬ 
ize  the  phase -only  null  synthesis  problem. 

We  conclude  the  report  by  demonstrating  an  equivalence  between  the  null 
synthesis  problem  treated  in  terms  of  obtaining  a  best  fit  to  a  given  set  of  element 
weights,  and  treated  in  terms  of  obtaining  a  best  approximation  to  tne  original  pattern. 
At  the  same  time  we  also  extend  the  null  synthesis  problems  studied  in  the  first  part 
of  the  report  by  considering  a  more  general  minimization  criterion  that  includes  the 
criteria  treated  earlier  as  special  cases. 


3.  ANALYSIS 

In  the  following  analysis  we  consider  a  linear  array  of  equispaced  isotropic 

elements  (see  Figure  1).  The  spacing  between  the  elements  is  d  and  the  phase 

reference  center  is  taken  to  be  the  center  of  the  array.  Lettirg  w  ,  n  ■  1, 2  ... ,  N, 

ttl  ** 

be  the  complex  weight  of  the  n  array  element,  the  array  field  pattern,  p(u)  is 
given  by 


(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  67.) 
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*  i  ■ 

N  N-«  M-2 


Figure  1.  Geometry  of 
Array 


p(u> 


N 
n*  1 


-Jdn« 


where 

dn  *  ■  •  (n  •  l)i  n  «  ta2|  >111  N 

and 


u  *  kd  sin  0 


with 


and  0  the  angle  measured  from  broadside  to  the  array.  *  Note  that  the  dfl  are  odd 
symmetric  with  respect  to  the  phase  reference  center;  that  is. 


V  -dN-n+r 


n  =  1, 2,  « •  • ,  N  • 


The  general  problem  we  address  here  is  the  following.  Suppose  vrt  are  given 
an  amplitude  taper  an>  n  •  1,2,  ....  N.  for  the  element  excitations  wr  example, 
for  low  sidelobea),  and  a  direction  ufl  for  the  peak  of  the  array  pattern  (look  direction). 
Then  the  array  coefficients  are 


Jd, 


on 


an  * 


nus 


n  =  1,  2,  . . . ,  N. 


The  amplitude  distribution  is  assumed  to  be  symmetrical  with  respect  to  the  array 
center.  We  now  wish  to  alter  the  element  weights  so  as  to  place  nulls  in  the  array 
pattern  at  a  set  of  M  prescribed  sidelobe  locations  u^,  k-  1,2  ....  M.  What  should 
the  coefficient  perturbations  be?  We  consider  two  principal  cases  of  this  problem: 
(1)  the  perturbations  are  of  both  amplitude  and  phase;  and  (2)  the  perturbations  are 
of  the  phases  only. 

*In  this  report,  we  use  the  implicit  time  dependence  e"^4**1. 
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2.1  Ajaplitud*  Mid  FHuc  Porturbation* 

We  begin  by  considering  the  case  in  which  both  the  amplitude  and  the  pnase  ol 
the  element  weights  can  be  perturbed  from  their  original  values.  The  perturbed 
coefficients  can  be  represented  as 


jd  u  jd  u 

w  »  a  e  n  8  +  an  e  n  8  (A  +ji  ). 
ft  n  n  n  n 


(u 


The  first  term  on  the  rigid  hand  side  of  Eq.  (1)  is  the  initial  value,  w  ,  of  the  n 
element  weight,  and  the  second  term  is  the  perturbation  of  the  weight 


th 


Aw  *  a  e  n  e 
n  n 


(An+iV- 


The  array  field  pattern  is  then 
N 

p(u> 


-  E  [*„  *  nUa  ♦  «n  *  +  i  ♦„>] 

N  .  . 

jd  u 

«  p  (u)  +  a  e  8  (A  +  j  0  >  e 
o  n  n  ”n 


-jdnu 


-jdnu 


(2) 


where  P0(u)  is  the  original  or  unperturbed  array  pattern.  We  now  wish  to  deter¬ 
mine  the  An  and  the  0n  to  place  nulls  at  the  M  locations  u^.  k«  1,  2,  . . . ,  M.  or 
equivalently  to  find  solutions  to  the  equation  system 

”  j u  u  jd  u. 

2^  an  (An  +  j  d>n)  c  “  8  e  n  «  -  P0<uk>,  k  »  1,  2 . M  .  (3) 


n=  1 


Note  that  in  this  equation  system  there  are  2N  unknowns  and  2M  equations 
(considering  the  real  and  imaginary  parts  of  Eq.  (d)  separately].  Hence,  if  the 
number,  M.  cf  prescribed  null  locations  is  less  than  N,  the  equation  system  as  It 
stands  does  not  have  a  uniquely  determined  solution.  Clearly  a  further  requirement 
(or  requirements)  must  be  imposed  on  the  solution  to  determine  it  uniquely. 

At  this  point,  motivated  by  the  desire  to  keep  the  weight  perturbations  as  small 
as  possible,  two  possibilities  suggest  themselves.  The  first  is  to  find  the  solution 
to  Eq.  (3>  which  minimizes  the  sum  of  the  squares  of  the  absolute  values  of  the 
weight  perturbations  relative  to  the  original  weights. 

E  |  A  +  J  A  | 2  -  E  (A  1+&. 

_  _  .  n  n  _ ,  n  rn 

n«=  1  n=  1 
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The  second  possibility  is  to  find  the  solution  which  minimizes  the  sum  of  the  squares 
of  the  absolute  values  of  the  total  weight  perturbations 


We  now  treat  each  of  these  two  cases  in  turn.  The  treatment  of  both  cases  is  based 
on  the  following  result  which  we  state  as  a  theorem: 

THEOREM  A: 

Let  A  be  an  m  X  n  complex  matrix  of  rank  m,  m  <  n  (that  is.  A  has  m  linearly 
independent  rows).  Then  the  solution  x  of  the  equation  system 

Ax  »  jr 

which  minimizes  the  norm  of  x.  |  |x|  |  *  (x^x)*^,  can  be 

expressed  as  a  linear  combination  of  the  columns  of  A^.  (Here  t  represents  the 
Hermitian  transpose  operator. ) 

PROOF: 

The  proof  is  an  immediate  corollary  of  the  fact  that  the  minimum  norm  solution 

is10 


X  -  At(AAt)‘1i. 

If  we  let  b  be  the  m  X  1  vector  (AA^)_1^then 

X  m  A^b 

so  that  x  is  a  linear  combination  of  the  columns  of  A^  with  coefficients  from  b. 
Applying  this  theorem  first  to  finding  the  solution  of  Eq.  (3)  which  minimizes 

<+  $  we  let  _ 

£-  J*t»  A2+  j*2.  ....  An  +  j*MJT  (4) 

X  P0(u2) . P0(um„T  (5) 


and 


10,  Rao.  C.,  and  Mitra,  S.  ( 197 1)  Generalized  Inverse  of  Matrices  and  Its 
Applications.  -  John  Wiley  &  Sons,  New  York,  p  46.  1 
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Jdl(u8-“l>  id2(uS-ul) 


...  aNe 


A*  a!« 


jdl(ua  *  U2>  jd2(u8‘u2>  aNe 


eJcVus  -  “l* 

JVus  ‘  V 


4  Jd l(us*V  n  id2<tt.-“M>  a  JdN<U8'uM>! 
al*  a2  e  aN*  I 


ao  that  the  system  of  Eqs.  (3)  can  be  written  as  Ax*  j.  It  follows  from  the  theorem 
that 


jd  <um  *  UE> 

A+j*=£baenm  8 
n  J  ya  m  n 

m  =  l 


.  n  =  1.2,...,  N 


where  the  vector  of  coefficients,  b,  is  given  by  inverting  the  equation 
(AA^)b  *  jr  . 

The  elements  of  the  MX  M  matrix  AA^  are  given  by 

t  Jdn(u«  *  V  *jd  (u  “  Um> 

<“ \m  .  £  %  *„« 


*  X)  cos  {d^u^  -  um)J»  k,  m  «  1, 2 . 

n  *  1 


making  usc  of  the  odd  symmetry  of  the  dft  and  the  assumed  even  symmetry  of  the 
original  distribution  of  the  element  weight  amplitudes.  The  elements  of  AA^  are 
thus  real.  The  components  of  the  vector  £,  the  negatives  of  the  unperturbed  pattern 
values  at  the  prescribed  null  locations,  are  likewise  real  since  the  unperturbed 
pattern  is  real  because  of  the  symmetry  of  the  dn  and  of  the  original  amplitude 
distribution.  Hence  the  coefficients  b^.  m  •  1, 2.  ...  M,  are  real  and  it  follows 
fiom  Eq.  (7)  that 


14 


M 

A  =  a  2_i  b  ccs  Id  (u  -  u  )] 
n  a  i  t  m  1  n  m  s' 

M 

0  =  a  ),  b  sin[d  (u  -  u  )]  . 

’  n  n  *—*.  m  1  n  m  s  1 

m  =  1 

The  A^  arc  thus  even  symmetric  and  the  $n  odd  symmetric  with  respect  to  a: 

that  is. 


^n  -  n+  1  j 

?n  =  ‘?N  n+  1  I 


n  =  1,2,  ....  N  . 


The  representation  of  the  relative  weight  perturbation  in  the  form  of  Eq.  (7) 
can  be  given  an  appealing  physical  interpretation.  Substituting  Eq.  (7)  into  Eq.  (2) 

we  obtain 


Aw  =a2  E  b  ejd"Um 

n  r  Jr?  \  rn 

m=  1 


and  thus 

p(u)  =  po<u)  +  Ap(u) 


I  ; 


where 


M 


Ap(u)  =  E  b 


N 

V 


2  -jdn(u 

-  .  a  e 
i  ni  n 

m  -  1  n  =  1 


U  ) 
m 


(9) 


Accordingly  the  perturbed  pattern  equals  the  original  pattern  plus  the  sum  of  M 
cancellation  beams  each  corresponding  to  a  taper  of  the  amplitude  of  the  element 
excitations  proportional  to  the  square  of  the  original  amplitude  taper  and  with  a 
pattern  peak  in  the  direction  of  one  of  the  desired  null  locations.  In  the  sequel  we 
will  sometimes  refer  to  the  cancellation  beams  corresponding  to  the  amplitude  taper 
an2  as  "tapered  cancellation  beams.  " 

N  ,  2  N  ,  ,»  N  -  -  - 

If  instead  of  minimizing  £  <A^  +  we  minimize  £  |Awn|  •  a*(A*+^*) 

(see  p.  13),  ttien  similarly  to  Eq.  (4)  and  Eq.  (6)  we  let 


i 
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*  “  [*j  j  +  J  *2  +  j  ^2^*  •  •  •  •  +  ^ 


and 


l'“a  “1'  J  “2'“»  “1'  JdN(u«'ul) 


jd  <u  -  u.)  J<2,(u«  -  u,> 

e  1  8  1  e 


jdx<ua-u2)  Jd2<ua  *  u2*  JV.'V 

e  e  . .  •  e 


Jdl(tt.-UM)  Jd2(ua  -  V  JdN<u.*uM) 
e  e  4  . .  e 


thereby  obtaining 


M 


V  jd(u  -  u  ) 

(A  +  j  0  ) 3  2-  be  n  m  n-1,3.  ....  N 


n"  n  *  rn  rr.  m 
m=  l 


The  vector  of  coefficients,  b,  ia  obtained  aa  before  by  inverting  Eq. 
elements  of  the  matrix  AA^  are  now  given  by 
N 

[aa*]  *  53  cos  Jdn  (i^  -  um)J.  k,  m  •  1,2,  ....  M  . 

The  coefficients  b  .  m  =  1, 2,  ....  M  are  again  real  and  we  obtain 

m 

M 


A  =  —  ^  b  cos  [d  (u  -  u  )] 

n  a_  m  n  m  s' 

n  m=  l 

M 

6  =  ~  J3  b  aln(d(u  -  u  )]  . 

rn  a  m  n  m  s' 

n  in  -  l 


The  weight  perturbation  takes  the  form 

t — '  j  d  u 

Aw  »  /  ,  b  e  ,  n  «  1, 2,  . . . ,  N 

n  ,  m 

m  =  1 


and  the  pattern  perturbation  ia 


(81  where  the 


U0> 
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(11a) 


Ap(u)  *  -5^  b  ]C 


-jd  (u  -  u  ) 
.  n  m 


m  =  1  4-1  n*  1 


sin(N(u  -  u  >/ 2] 


m  sin  Hu  -  u  J 
m  *  i  1  m 


The  perturbed  pattern  is  thus  equal  to  the  sum  of  the  original  pattern  and  M  cancella- 
tion  beams  each  corresponding  to  a  uniform  amplitude  distribution  of  the  element 
excitations  and  with  a  peak  in  the  direction  of  one  of  the  null  locations.  We  will 
sometimes  refer  to  such  cancellation  beams  as  "uniform  cancellation  beamp  ' 

The  representation  of  the  total  pattern  perturbation  as  the  sum  of  M  cancella¬ 
tion  beams  corresponding  to  uniform  element  amplitudes  is  alro  the  form  of  the 

pattern  perturbations  for  a  half -wave  length -spacing  linear  a  rr  »y  obtained  by 
o 

Steyskal  who  finds  the  weights  that  will  give  the  pattern  with  nulls  at  a  set  of  pre¬ 
scribed  locations  which  differs  as  little  as  possible  in  the  mean  square  sense  from 
the  original  pattern.  In  Appendix  A  we  give  an  alternate  dei  ivetion  of  this  result 
as  a  special  case  of  the  solution  to  the  minimum  pattern  perturbation  problem  for  a 
linear  array  with  arbitrary  uniform  element  spacing.  Thus  the  solution  to  the  prob¬ 
lem  of  placing  nulls  at  a  prescribed  set  of  locations  which  minimizes  the  aim  ot 

the  squares  of  the  absolute  values  of  the  total  weight  perturbations 
N  -  N  ,  .  , 

TJ  I  Aw  I  »  V  a  (A  +  $  )  is  also  the  solution  that  minimizes  the  mean 
n*=  I  n  n=lnn  n 

square  difference  between  the  original  field  pattern  and  the  perturbed  field  pattern. 

In  low  sidelobe  applications  we  can  anticipate  that  the  minimum  (in  either  of  the 
two  cases  we  have  treated  above)  weight  perturbations  required  to  place  nulls  at  a 
set  of  sidelobe  locations  will  be  small.  Starting  with  the  expression  for  the  ntJl 
perturbed  weight 


*  *  •_  e 

n  n 


(1  +  An+J*n> 


1  d  u_  2  2 

*n  e  l(1+V  +*n> 


1/2  itwn'i(-LTzr) 


we  can  then  expand 


[‘*i  (i%rf  ♦•••] 
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and  retain  only  the  linear  terms  thus  obtaining 


id  u  id 

*  (1  ♦  A  )  e  n  *  e  V 
.n  n 


Hence  (or  small  A_  a  id  6  ,  A  Is  the  relative  change  In  the  amplitude  and  6  is  the 
n  nth  n  n 

phase  perturbation  of  the  n  weight. 


2.2  Phan-Only  Perturbation* 


We  now  turn  our  attention  to  the  case  in  which  phase  perturbations  only  are 
allowed  in  the  element  weights.  The  amplitudes  of  the  element  excitations  are  fixed 
at  their  original  values.  The  phase-only  analysis  parallels  that  performed  above  in 
which  both  amplitude  and  phase  were  allowed  to  vary  from  their  initial  values. 

The  perturbed  coefficients  are  represented  by 


jd  u 

an*  “  *  * 


(lla) 


Jdnus 

ane  +  an 


jd  u  Id 
e  n  8  <eJ9a 


1). 


1.2. 


N 


(12b) 


We  assume  that  the  phase  perturbations  are  small  (for  example,  for  applications  to 
null  synthesis  in  low  sidelobe  array  patterns).  Then 


1  "j^n 


(13) 


and 


j u  jd  u 

_  M  .  _  is  •  .  j  _  a  n  8 

w  *  a  e  +  J  *.  e 
n  n  n  n 


(14) 


The  array  pattern  is  giver,  by 

«  Jd  u  -jd  u 

P(u)  *  PQ(u)  +  jZ#“n^nenSe  n 
n*  1 


and  the  equation  system  for  the  nulls  is 
N 

an  ain  'dn{uk  *  “a**  “  ”  Po^ 
N 

n?i  #n*nC08,dn<Uk-us)I  “  0 


k  -  1.2,  ....  M 


(15a) 


<15b) 
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This  equation  system  has  2  M  equations  for  N  unknowns  and  hence,  as  in  the  case  of 
perturbations  of  both  amplitude  and  phase,  a  unique  solution  to  tlie  system  of  equa¬ 
tions  can  be  obtained  only  if  a  further  requirement  is  imposed  on  the  solution.  (It 
is  assumed  of  course  that  the  number  of  unknown  phase  perturbations  is  greater 
than  twice  the  number  of  null  locations. )  As  above  we  treat  two  cases,  the  first 
here  corresponding  to  the  requirement  that  the  solution  minimise  the  sum  of  the 

N  2 

squares  of  the  phase  perturbations,  £  p*  (or  equivalently,  minimise  the  sum  of 

n*  1  n 

the  absolute  values  of  the  weight  perturbations  relative  to  the  original  weights),  and 

the  second  to  the  requirement  that  the  solution  minimise  the  sum  of  the  squares 

N  , 

of  the  absolute  values  of  the  total  weight  perturbations,  T'  (a  p  )  . 

n*  1  n  n 

Consider  first  minimising  the  sum  of  the  squares  of  the  phase  perturbations. 

We  write  the  system  of  Eqa.  (IS)  in  the  form 

| 

Ax  -  jr  i  (18) 


where 


1“  *2*  •••» 


‘sj  sin(d1<u1  -  uB»  aa  sintdj^j  -  ug)J  . . .  Sg  sinld^Uj  -  u^)] 

! 

i 

at  statdj  (u2  -  us)]  a2  sln(d2(u2  -  us)J  ...  s^  sin(dN(u2  -  ug)J 


*1  •l“lVuM  •  us’l  *2  ,la<d2(uM  ■  VI  —  *N*toldN(uM  -  us,)| 


and  jr  is  given  by  Eq.  (5).  Using  Theorem  A  we  then  obtain  the  representation  of  the 
phase  perturbation  as 
M 

♦n"*n  £  bm  ■ia[dn<um  “  V1  * 

m*  1 

The  vector  of  coefficients,  b,  is  obtained  by  inverting  the  equation 


(AAa )b  ■  jr 
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with  the  elements  of  the  matrix  AA1  given  by 

[^Ikm  "  n?i  ■»  ‘inIdn  <Uk  ‘  Ua>1  ,lnldn<Um  *  Us 


(21) 


It  follows  from  Eq.  (19),  because  of  the  even  symmetry  of  the  amplitudes  an 
assumed  throughout  and  the  odd  symmetry  of  the  dn»  that  the  are  odd  symmetric 
with  respect  to  the  phase  reference  center;  that  Is, 

*n*-*N-n+r  n“1*2 . N- 


This  means  that  Eq.  (ISb)  is  automatically  satisfied  by  the  solution  obtained  to 
Eq.  (15a). 

It  is  important  to  note  that  the  values  of  pn  calculated  from  Eq.  (19)  must  be 

substituted  in  Eq.  (12a)  or  Eq.  (12b)  to  obtain  the  perturbed  weights.  Since  the 

values  for  the  0n  were  obtained  using  the  approximation  (14),  for  the  weights, 

this  means  that  the  resulting  pattern  nulls  cannot  be  expected  to  be  theoretically 

perfect  as  they  would  be  if  the  calculated  values  of  a  were  substituted  in  approxi- 
♦  ** 

mation  (14). 

The  approximated  form  (14)  la  useful,  however,  for  obtaining  a  physical  inter¬ 
pretation  of  the  phase-only  perturbation  solution.  Substituting  Eq.  (19)  in  (14)  we 
obtain  the  approximate  weight  perturbation 


Aw 


w  -  a  e 
n  n 


Jdnus 


i*n 

^n  * 

M 

J»n 

E 

m»  1 

„2 

M 

an 

T 

E 

m=  1 

Jdnus 


J  dn  us 


bme  “  •«ln(dn(um-u>)] 


The  value  of  the  perturbed  pattern  at  a  prescribed  null  location  u^  is  in  fact  given  by 

y  .  IJ+n  ,.  .  iA  ,1  ~^dn “k  y  *n  "^“k 

n- t  »  L*  '  (l  +  i  ♦„»  J  e  *  n~  an  T  e 

with  the  value  of  the  power  pattern  at  the  location  given  approximately  by 

1  y  y  .  _  .2  2  *J(dn'dm)uk 

*n£l 
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ar.d  the  corresponding  pattern  perturbation 

V'  an  T  “j®L(u  -  u  )  -jd  (u  -  {  2u  -  u  \  )"| 

*ip(u)  *  2]  bm  S  -7-  J  *  m  -e  “  '/J  <22) 

m  -  1  n=  1  L 

The  perturbed  pattern  is  thus  approximately  equal  to  the  sum  of  the  original  pattern 
and  the  sum  of  the  differences  of  M  pairs  of  cancellation  beams.  Each  of  these  pairs 
of  beams  corresponds  to  an  element  amplitude  taper  proportional  to  the  square  of 
the  original  amplitude  taper,  with  one  member  of  the  pair  having  a  peak  in  the 
direction  of  a  null  location  and  the  other  member  of  the  beam  pair  having  a  peak  in 
the  symmetric  direction  on  the  other  side  of  the  original  pattern  look  direction,  ug. 
[Note  that  2u  •«  ■  u  -  2  (u_  -  u  ).  I  Titus  the  placing  of  nulls  in  the  prescribed 

locations  is  accompanied  by  a  doubling  of  the  pattern  at  symmetric  locations  on  the 
other  side  of  the  look  direction. 

The  requirement  of  minimizing  the  sum  of  the  squares  of  the  absolute  values  of 

N  , 

the  total  weight  perturbations  £  (*„♦„>  »  (•“  P-  19)  leads,  similarly  to  F.q.  (17) 

n*  1  n  n 

and  Eq.  (18).  first  to  expressing  the  system  of  Eqs.  (15a)  in  the  form  of  Eq.  (16) 
where  now 


x  -  [aj0r  »2 #2* 


■sinfdjfaj  -  ug)l  sinldjtaj  -  ug)[  ...  sinIdjj'Uj  -  ug)l‘ 

sln(dj  (Uj  -  ug)J  sln[d2(u2  -  u#)l  . . .  sin  [d^  (t^  -  ug)[ 


sln[dj (Um  •  ug)|  sin[d2(uM  -  ug)J  ...  •lnldN<uM*V) 


Hence,  similarly  to  Eq.  (19;, 


m 

■  ♦_  •  53  *>_  *in[d  (u  -  U,M 

n  n  m  n  m  s 


with  the  vector  of  coefficients,  b,  obtained  by  inverting  Eq.  (20)  with 


21 


N 

E,  .lnldn^  -«.)]  •i«»[dn(ttm  -u,)]  . 

f\  ■  X 


The  weight  perturbation  is  then  approximated  by 

.  “  f  j d  u  jd  <2u  -u  >1 

b»L*  -  J 


and  the  corresponding  pattern  perturbation  by 


•  V  bm  [e*Jdn(u-um)  {2u.‘um()] 

m  •  1  "S’  L  *  J 


-  E 

m  *  1 


sin(N(u-um)/2)  sin[N(u  -  |2u(-um  |)/2) 

sin[(u  -u  >/ 2J  <2u.  -u  I  '  -i 

m  sin  (u  -  |  *  m  f  )/2  I 


.  (24b) 


Minimising  the  squares  of  the  total  weight  perturbations  thus  leads  to  a  result 
similar  to  that  for  minimising  the  squares  of  the  phase  perturbations  (relative 
weight  perturbations).  The  perturbing  or  sidelobe  cancelling  pattern  is  again 
approximated  by  the  sum  of  the  differences  of  M  pairs  of  beams,  each  of  which 
corresponds  here  to  a  unifo-m  amplitude  distribution  of  the  element  weights.  Equa¬ 
tion  (23)  shows  that  the  phase  perturbations  are  inversely  proportional  to  the  ampli¬ 
tudes  of  the  weights,  whereas  for  the  minimum -phase -perturbation  solution. 

Eq.  (19).  the  phase  perturbations  are  directly  proportional  to  the  weight  amplitudes. 
It  is  also  simple  to  show  (see  Appendix  A)  that  for  half -wave length  element  spacing, 
the  pattern  given  by  the  approximation  (24)  which  corresponds  to  minimizing 
the  approximate  total  weight  perturbations,  is  also  the  pattern  that  differs  least  from 
the  original  pattern  in  a  mean  square  sense,  while  having  nullB  at  the  desired  loca¬ 
tions  obtained  with  small  phase-only  perturbations. 

It  is  instructive  to  compare  the  approximate  (for  small  phase-only  perturbations) 
representations  of  the  pattern  perturbations  obtained  here,  (22)  and  (24),  with 
the  representations  in  Eqs.  (9)  and  (11),  obtained  earlier  when  both  the  amplitude 
and  phase  of  the  weights  are  allowed  to  vary.  Minimizing  the  sum  of  the  squares  of 
thj%ab"°lUte  v*lues  ot  th*  relative  (to  the  original  amplitude  weight)  perturbations 

<*  n  -  I  »  jfn  for  phase-only  perturbations  and  An  +  j+n  for  combined  amplitude 
and  phase  perturbations)  leads  to  a  representation  of  the  cancellation  pattern  in 
terms  of  beams  corresponding  to  an  element  amplitude  taper  proportional  to  the 
square  of  the  original  amplitudes;  whereas  minimizing  the  sum  of  the  squares  of 

Ju  ju 

the  absolute  values  of  the  total  weight  perturbations  (a  e  “(e  n-l)*jae  * 
j  us  n  *  n  * 

*„  *  (£„  ♦  respectively)  leads  to  a  representation  of  the  cancellation  pattern 
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in  term*  of  beams  corresponding  to  a  uniform  amplitude  distribution.  Phase-only 
perturbat'ons  require  taro  beams  for  each  desired  null,  one  pointing  in  the  direction 
of  the  null  and  the  other,  of  opposite  sign,  in  the  symmetric  direction  on  the  opposite 
side  of  the  original  main  beam  axis,  with  the  result  that  the  placement  of  the  nulls  is 
achieved  at  the  expense  of  doubling  the  sidelobe  levels  at  the  symmetric  locations. 
When  both  amplitude  and  phase  are  allowed  to  vary,  there  is  only  one  beam  for  each 
null  pointing  in  the  direction  of  the  null  location.  For  half-wavelength  element 
spacing,  the  patterns  corresponding  to  minimizing  the  sum  of  the  squares  of  the 
total  weight  perturbations  are  also  the  patterns  that  differ  least  in  the  mem  square 
sense  from  the  original  patterns.  Because  of  the  small  phase  approximation  (13) 
employed  in  the  phase-only  analysis,  the  resulting  pattern  nulls  are  not  perfect, 
whereas  the  patterns  achieved  with  both  amplitude  and  phase  control  have  theoretic¬ 
ally  perfect  nulls  at  the  desired  locations. 

Before  discussing  numerical  results,  it  is  worthwhile  to  introduce  a  useful 
normalization  of  the  beam  coefficients.  Consider  first  Eqs.  (9)  and  (11)  obtained 
for  the  perturbing  or  sidelobe  cancelling  pattern  Ap(u)  when  both  the  amplitude  and 
the  phase  of  the  weights  are  allowed  to  vary  from  their  original  values.  For  simpli¬ 
city  suppose  there  is  only  one  null  location  (that  is.  1).  Then  Eq.  (9)  gives  for 
minimum  squared  relative  weight  perturbations, 

N 

-Poiu^-  Ap(Uj)  ■  bj  2 
n*  1 


so  that 
b. 


Po(ul) 

~~N 

z  ■: 

n«  1  n 


whereas  Eq.  (11)  gives,  for  minimum  squared  total  weight  perturbations, 
■P0<u.>  *  AP(“j>  *  Nbj  \ 


or 

b  -  P°<Ul> 

bl  *  N - 

The  beam  coefficients  in  the  two  cases  can  have  widely  differing  values  depending 
on  the  amplitude  taper,  an.  and  the  number  of  elements  in  the  array,  NJ  If  simul¬ 
taneous  nulls  at  multiple  locations  are  desired,  then  the  deviations  of  the  values  of 
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the  beam  coefficients  from  tneir  values  obtained  for  the  aame  null  locations  taken 
one  at  a  time  give  a  measure  of  the  extent  of  coupling  or  interference  between  the 
component  beams  of  the  perturbing  pattern.  For  ease  in  describing  the  extent  of 
coupling  between  component  beams  as  well  as  for  comparing  differences  in  coupling 
corresponding  to  different  types  of  beams ,  it  is  desirable  to  normalize  the  beam 
coefficients  so  that  their  value  equals  the  negative  of  the  unperturbed  pattern  value 
at  the  null  location  when  there  is  only  one  null  location.  This  is  equivalent  to  re¬ 
quiring  that  each  component  beam  of  the  cancelling  pattern  have  a  magnitude  of 
un».y  in  the  peak  of  direction  of  the  beam.  It  follows  *hat  the  beam  coefficients  for 

"tapered  cancellation  beams"  she  ild  be  multiplied  by  ^  a2,  and  the  beam  coeffi- 

n*  1  n 

cients  for  "uniform  cancellation  beams"  should  be  multiplied  hy  N.  The  two  cases 
can  be  given  a  uniform  representation  as  follows.  Let 


a^,  tapered  cancellation  beams 


1,  uniform  cancellation  beams 


Then  the  equation  system  for  the  normalized  beam  coefficients,  b^,  is  (see  Eq.  (8)] 
(AA^)  b  *  £ 


where 


w 

[aA^J  *  53  tn  cos{dn(uk  -  k,  m  ■  1, 2,  ....  M  . 


km  n=  1 


The  n  weight  perturbation  is  given  by 
A.  .<  Zb' 

n  n  m  \ 


and  the  cancellation  pattern  by 
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V  '  ^  •■>dn<u-um) 

AP(u)  =  E  bm  E  ta 


The  real  and  imaginary  parts  of  the  relative  weight  perturbations  are  given  re¬ 
spectively  by 


4 

M 

n 

a 

E 

bm 

cos (dn 

(Um  *  VI 

m  s 

n 

m=  1 

t 

M 

n 

au 

E 

m=  1 

t 

b 

m 

sin  [d 

1  n 

(u  - 
m 

The  normalization  for  phase-only  control  is  entirely  similar.  With  tn  defined 
by  Eq.  (26).  the  equation  system  for  the  normalized  beam  coefficients  is  [cf.  Eq.  (20)] 

<AAT)  b'  *  £ 

where 


r  l  N 

Jkm  *  Z*  2,n  Bin<dn(uk  *  u.»  slnldn(um  *  U»)J  * 


The  phase  perturbations  are 
.  M 


4  »  2  -2-  E  b'  sin  Id  (u  -  u  ))  . 
n  *n  m*i  m  n  m  8 


The  perturbation  of  the  n  weight  can  be  represented  approximately  by 


Aw_  « 
n 


.  E  b' 

n  m=  1  m  L  J 


and  the  cancellation  pattern  approximated  by 


**.).  E  bm  E  ’l. 

tail  m  n«l  n  L  J 


1  NUMERICAL  RESULTS  AND  DISCUSSIONS 


In  this  section  we  present  and  discuss  the  results  of  computations  performed  to 
display  the  basic  features  of  the  solutions  obtained  in  the  previous  section.  All 
computations  were  dons  for  an  array  with  41  isotropic  elements  with  half  wavelength 
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spacing.  Unless  otherwise  stated,  the  original  or  unperturbed  pattern  corresponds 
to  a  40  dB  Chebyshev  taper  of  the  element  excitations.  *  * 

We  begin  by  discussing  the  simplest  case— that  of  imposing  a  single  null  in  the 
pattern.  A  location  of  8  =  15.  23*  was  chosen  corresponding  to  the  approximate 
location  of  the  peak  of  the  fifth  sidelobe  of  the  unperturbed  pattern.  Figures  2a-2c 
and  3a-3c  show  the  patterns  corresponding  to  perturbations  of  both  the  amplitudes 

and  phases  of  the  element  weights  minimizing  the  relative  weight  perturbations 
N  ,  N  „ 

V(|aw  |/a  )  [Eq.  (9)]  and  the  absolute  weight  perturbations,  £  lAw  I  , , 
n=  1  n  n  n=  1  . 

(Eq.  (11)J,  respectively.  For  each  of  these  two  cases  ,.e  show  three  graphs,  the 
first  with  the  perturbed  pattern  for  6  from  -90*  to  90",  the  second  with  the  can¬ 
cellation  pattern  from  -90®  to  90“.  and  the  third  with  both  the  original  and  the  per¬ 
turbed  patterns  in  the  range  from  0“  to  25  “. 

We  see  that  the  cancellation  beam  in  Figure  2b  for  minimized  relative  weight 
perturbations,  which  corresponds  to  a  taper  of  the  element  excitations  equal  to 
the  square  of  the  original  amplitude  taper,  has  a  much  broader  mainlobe  and  much 
lower  sidelobes  than  the  cancellation  beam  in  Figure  3b  which  corresponds  to  a 
uniform  distribution  of  the  element  amplitudes.  This  means  that  the  original  p?tt'  ,-n 
is  distorted  more  strongly  in  the  vicinity  of  the  null  for  minimized  relative  weight 
perturbations  than  it  is  for  minimized  total  weight  perturbations,  and  conversely 
that  the  original  pattern  is  distorted  less  strongly  far  away  from  the  null  for  the 
minimized  relative  weight  perturbations  than  it  is  for  minimized  total  weight  per  ¬ 
turbations.  (In  a  mean  square  sense,  as  is  demonstrated  in  Appendix  A,  the  dis 
tortion  of  the  original  pattern  is  minimized  by  the  minimized  total  weight  perturba¬ 
tion  scheme. )  Note  that  the  broad  mainlobe  of  the  cancellation  pattern  in  Figure  2b, 
being  out  of  phase  with  the  fifth  sidelobe  of  the  original  pattern,  la  in  phase  with  the 
fourth  and  sixth  sidelobes  of  the  original  pattern  with  the  result  that  these  sidelobe 
levels  are  raised  (see  Figures  2a,  2c).  Indeed  the  effect  of  the  mainlobe  of  this 
cancellation  pattern  extends  even  to  the  third  and  seventh  sidelobes  of  the  original 
pattern  which,  being  out  of  phase  with  the  mainlobe  of  the  cancellation  pattern,  are 
lowered  slightly  from  their  original  values.  By  the  time  the  second  or  eighth  side¬ 
lobes  are  reached,  however,  the  effect  of  the  cancellation  pattern  of  Figure  2b  is  no 
longer  felt  and  the  perturbed  and  original  patterns  are  virtually  indistinguishable. 

In  contrast,  the  cancellation  beam  of  Figure  3b,  which  tapers  much  more  slowly  than 
does  the  cancellation  beam  of  Figure  2b,  can  be  seen  by  referring  to  Figure  3c  to 
have  a  small  effect  as  far  as  the  first  sidelobe  of  the  original  pattern. 


11.  Drane,  C. ,  Jr.  (1964)  Dolph-Chebyshev  excitation  coefficient  approximation 
IEEE  Trans.  Antennas  Propae. ,  AP-1 2:781-782. 
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POWER  (OB) 


THETA  (DEGREES) 


Figure  2a.  Perturbed  Pattern  With 
One  Null  Imposed  at  15.  23*  With 
Combined  Amplitude  and  Phase 
Perturbations,  Minimizing  Relative 
Weight  Perturbations.  0  ■  >90*  to  +  90* 


Figure  2b.  Cancellation  Patterns  to 
Impose  One  Null  at  15.  23*  With  Combined 
Amplitude  and  Phase  Perturbations, 
Minimizing  Relative  Weight  Perturbations. 
0  *  -90*  to  +90* 


THETA  (DEGREES) 


Figure  2c.  Original  Pattern  and 
Perturbed  Pattern  With  One  Null 
Imposed  at  15.  23*  With  Combined 
Amplitude  and  Phase  Perturbations, 
Minimizing  Relative  Weight  Pertuba- 
tions.  0  -  0*  to  25* 
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THETA  (OECREES) 


Figure  3a.  Perturbed  Pattern  With 
One  Null  Imposed  at  IS.  23*  With 
Combined  Amplitude  and  Phase 
Perturbations,  Minimizing  Total  Weight 
Perturbations.  8  «  ■  90*  to  +90* 


Figure  3b.  Cancellation  Pattern  to  Impose 
One  Null  at  15.  23*  With  Combined 
Amplitude  and  Phase  Perturbations, 
Minimizing  Total  Weight  Perturbations. 

8  -  -90*  to  +90* 


THETA  (DEGREES) 


Figure  3c.  Original  Pattern  and 
Perturbed  Pattern  With  One  Null 
Imposed  at  15.  23*  With  Combined 
Amplitude  and  Phase  Perturbations. 
Minimizing  Total  Weight  Perturba¬ 
tions.  $  =  0*  to  25* 
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Figures  4a-4c  and  5a-5c  show  patterns  for  a  single  imposed  null  at  15.  23” 
corresponding  to  perturbations  of  only  the  phases  of  the  element  weights.  Fig¬ 
ures  4a-4c  show  the  results  when  the  sum  of  the  squares  of  the  phase  perturbations. 

N  2 

7*  is  minimized.  (Eq.  (22)],  and  Figures  5a-5c  chow  the  results  when 
n=  l 

N  2 

£  (a  $  ;  is  minimized  [Eq.  (24)].  The  principal  change  in  these  patterns  from 
n=  1 

those  in  Figures  2a-2c  and  3a-3c  corresponding  to  perturbations  of  both  the  ampli¬ 
tudes  and  the  phases  of  the  weights  is  that  now  the  cancellation  pattern  is  the  super¬ 
position  of  two  beams,  one  with  a  peak  in  the  direction  of  the  null,  and  the  other  with 
a  peak  of  opposite  phase  in  the  symmetric  direction  on  the  other  side  of  broadside. 
Hence  any  change  in  the  pattern  on  the  null  side  of  broadside  (8  >  0)  is  accompanied 
by  an  equal  and  opposite  change  on  the  other  side  of  broadside  (9  <  0).  Note  es¬ 
pecially  that  the  sidelobe  peak  at  -15.  23“  in  Figures  4a  and  5a  is  raised  S  dB.  This 
is  because  the  cancellation  pattern  to  produce  a  null  at  +15.  23“  must  have  a  magni¬ 
tude  equal  to,  and  a  phase  opposite,  that  of  the  original  pattern  at  15.  23“  and  hence 
adds  in  phase  to  the  original  pattern  (which  is  symmetric  around  8  =  0)  at  -15.  23“ 
doubling  its  magnitude  there. 

The  cancellation  pattern  of  Figure  4b  is  the  superposition  of  two  beams  of  the 
form  of  Figure  2b  (that  is,  beams  corresponding  to  an  element  amplitude  taper 
equal  to  the  square  of  the  original  taper),  while  the  cancellation  pattern  of  Figure  5b 
is  the  superposition  of  two  beams  of  the  form  of  Figure  3b  (that  is,  beams  corres¬ 
ponding  to  a  uniform  element  amplitude  distribution).  On  the  null  side  of  broadside, 

0  >  0,  the  phase-only  cancellation  patterns  are  dominated  by  exactly  the  same  beams 
that  are  used  for  cancellation  in  the  corresponding  patterns  when  both  amplitudes 
and  phases  of  the  element  excitations  vary.  Hence  the  basic  features  of  the  phase- 
only  patterns  for  positive  8  are  the  same  as  those  noted  above  for  combined  ampli¬ 
tude  and  phase  perturbations.  Figures  4a-4c  being  similar  to  Figures  2a-2e  and 
Figures  5a-5c  being  similar  to  Figures  3a-3c. 

In  Table  1  we  have  tabulated  the  values  of  the  beam  coefficients,  and  the  values 
of  the  original  and  perturbed  patterns  at  the  null  positions  for  this  and  succeeding 
cases.  For  one  nuU  at  15.  23“  and  perturbations  of  both  amplitudes  and  phase,  it 
wUl  be  noted  that  the  value  of  the  beam  coefficient  equals  the  negative  of  the  original 
pattern  at  the  null  position.  This  is  a  consequence  of  the  normalization  of  the  can¬ 
celling  beams  discussed  at  the  end  of  the  previous  section.  For  nhase-only  perturba¬ 
tions,  the  beam  coefficient  correspending  to  the  minimization  of  the  sum  of  the 


squares  of  the  phase  perturbations  aga.'n  equals  the  negative  of  the  original  pattern, 

N 

but  the  beam  coefficient  correspondtig  to  the  minimization  of  £ 


(an^n>  is  slightly 


larger.  This  deviation  of  the  latter  beam  coefficient  is  a  consequence  of  the  fact 
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that  for  phase-only  perturbations  and  one  imposed  null,  che  beam  directed  at  the 
null  position  must  not  only  cancel  the  original  pattern  there  but  must  also  cancel 
the  sidelobe  of  the  beam  directed  at  the  symmetrical  position  on  the  other  side 
of  broadside.  Since  the  sidelobes  of  the  beams  used  for  cancellation  when 
N  2 

Y.  0  is  minimized  taper  off  considerably  more  quickly  than  the  sidelobes  of  the  beams 
11=1  “  N  2 

used  for  cancellation  when  Y  <a  0  '  *a  minimized  (compare  Figures  2b  with 

n=  1  n  n 

Figure  3b),  the  effect  of  the  sidelobe  of  the  beam  directed  at  the  symmetrical  posi¬ 
tion  of  the  null  is  more  strongly  felt  In  the  latter  case  than  it  is  in  the  former. 

Note  also  that  the  depth  of  the  nulls  achieved  when  both  the  amplitudes  end 
phases  of  the  element  coefficients  are  allowed  to  vary  are  considerably  deeper  than 
the  null  depths  when  only  the  phases  are  allowed  to  vary.  This  is  a  consequence  of 
the  small  angle  approximation  (13)  discussed  in  the  previous  section,  used  to  linear¬ 
ize  the  phase-only  nulling  problem.  The  phase  perturbations  obtained  by  solving  the 
linearized  nulling  problem  do  not  give  exact  nulls  when  substituted  back  into  the 
element  weights.  However,  the  smaller  the  phase  perturbations,  the  better  the 
approximation  is.  This  explains  why  the  null  depth  achieved  when  the  sum  of  the 

squares  of  the  phase  perturbations  is  minimized  is  deeper  than  the  null  obtained 

N  2 

when  the  weighted  sum  Y  (a  )  is  minimized— the  phase  perturbations  in  the 
n=  1 

former  case  are  in  a  least  squares  sense  smaller  than  those  in  the  latter  case. 

Following  this  discussion  of  imposing  a  single  null,  we  now  look  at  the  case  of 
two  null3  imposed  symmetrically  at  the  -3  dB  points  of  a  sidelobe.  Locations  of 
8  =  14.  54*  and  15.  94“  were  chosen  corresponding  to  the  -3  dB  points  of  the  fifth 
sidelobe  of  the  unperturbed  pattern.  Figures  6a-6c  and  7a-7c  show  the  pattern 
corresponds  to  the  perturbations  of  both  the  amplitudes  and  phases  of  the  element 
excitations  that  minimize  the  sum  of  the  squares  of  the  relative  weight  perturbations 
and  of  the  total  weight  perturbations  respectively.  Since  the  points  at  which  the  nulls 
are  imposed  belong  to  the  same  sidelobe  and  are  approximately  symmetrical  with 
respect  to  the  sidelobe  maximum,  the  beam  coefficients  are  approximately  equal 
(see  Table  1).  The  beams  appear  to  add  constructively  over  most  of  their  range  to 
form  what  has  very  much  the  appearance  of  a  single  beam  (compare  Figures  6b  and 
7b  to  Figures  2b  and  3b  respectively).  The  qualitative  features  of  the  patterns 
shown  in  Figures  6a-6c  and  7a-7c  are  entirely  similar  to  those  of  Figures  2a-2c 
and  3a-3c  respectively. 
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Figure  4a.  Perturbed  Pattern  With 
One  Null  Imposed  at  IS.  23*  With 
Phase-Only  Perturbations  Minimizing 

Xd  2.  0  ■  -90*  to  +90* 
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Figure  4b.  Canceilation  Pattern  to  m 

Impose  One  Null  at  15.  23*  With  Phase-Only  * 

Perturbations.  Minimizing  JO  2.  ^ 

9  -  -90*  to  +90*  n 
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Figure  4c.  Original  Pattern  and 
Perturbed  Pattern  With  One  Null 
Imposed  at  15. 23*  With  Phase-Only 

Perturbations,  Minimizing  T6  2. 

9  m  0*  to  25*  n 
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Figure  6b.  Cancellation  Pattern  to 
Impose  Two  T  ills  at  14.54*  and  15.94* 
With  Combined  Amplitude  and  Phase 
Perturbations.  Minimizing  Relative  Weight 
Perturbations.  0  *  -90*  to  +90* 
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Figure  6a.  Perturbed  Pattern  With 
Two  Nulls  Imposed  at  14.54*  and 
15.  94*  With  Combined  Amplitude  and 
Phase  Perturbations.  Minimizing 
Relative  Weight  Perturbations. 

6  =  -90*  to  +90* 


Figure  6c.  Original  Pattern  and 
Perturbed  Pattern  With  Two  Nulls 
Imposed  at  14.54*  and  15.94*  With 
Combined  Amplitude  and  Phase 
Perturbations.  Minimising  Relative 
Weight  Perturbations.  0  «  0*  to  25* 


j 


i 

Figure  7a.  Perturbed  Pattern  With  j 

Two  Nulls  Imposed  at  14.54*  and  ! 

15.  94*  With  Combined  Amplitude  and 
Phase  Perturbations,  Minimizing  : 

Total  Weight  Perturbations.  ! 

9  n  -80*  to  +90*  ) 


Figure  7c.  Original  Pattern  and 
Perturbed  Pattern  With  Two  Nulls 
Imposed  at  14.54*  and  15.94*  With 
Combined  Amplitude  and  Phase 
Perturbations,  Minimizing  Total 
Weight  Perturbations.  9  *  0*  to  25* 
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Figures  8a-8c  and  9a-9c  show  the  patterns  obtained  for  the  cases  of  the  same 

two  null  locations  at  the  >3  dB  points  of  the  fifth  sidelobe  when  only  perturbations 

of  the  element  phases  are  allowed.  Figures  8a-8c  correspond  to  the  min'mizatfon 

of  the  sum  of  the  squares  of  the  phase  perturbations,  and  Figures  8a-9c  to  the 
N  , 

minimisation  of  £  (a  *  ).  Referring  to  the  cancellation  patterns  shown  in 
n»  1  n  “ 

Figures  8b  and  9b  and  comparing  these  with  the  cancellation  patterns  in  Figures  4b 
and  5b  respectively,  we  see  that  the  two  pairs  of  beams  composing  the  phase -only 
cancellation  patterns  appear  to  add  constructively  to  form  one  resultant  beam  pair. 
The  qualitative  features  of  the  patterns  shown  in  Figures  8a-8c  and  9a-9c  are  thus 
the  same  as  those  noted  above  in  the  case  of  only  one  constrained  null  for  Figures 
4a-4c  and  5a-5c  respectively. 

Referring  to  Table  1  we  note  that,  unlike  the  case  of  one  null,  the  beam  coeffi¬ 
cients  for  all  four  minimization  types  are  considerably  smaller  than  the  corres¬ 
ponding  values  of  the  magnitude  of  the  original  pattern  at  the  null  locations.  This 
of  course  is  the  result  of  the  fact  that  the  cancellation  of  the  original  pattern  at  the 
null  locations  is  effected  by  the  peak  of  the  cancelling  beam  directed  at  one  null 
position  along  with  the  mainlobe  of  the  cancelling  beam  directed  at  the  other  null 
position  so  that  the  two  beams  help  one  another.  This  also  explains  why  the  beam 

N  2 

coefficients  obtained  when  £  (I  Aw  I  /a)  is  minimized  are  smaller  than  those 
N  n* 1  n  « 

obtained  when  Y\  I  Aw  I  2  is  minimized.  The  mainlobe  of  the  cancellation  beams 
n- 1  n 

in  the  former  case  falls  off  considerably  more  slowly  than  does  the  mainlobe  in  the 
latter  case  (compare  Figure  2b  with  Figure  3b)  and  so  the  peak  of  the  beam  directed 
at  one  null  position  can  get  more  assistance  in  cancelling  the  original  pattern  from 
the  mainlobe  of  the  beam  directed  at  the  other  null  position.  Some  of  the  same 
features  noted  above  for  the  case  of  one  null  are  seen  here  as  well.  Combined  ampll 
tude  and  phase  perturbations  give  deeper  nulls  than  do  phase  perturbations  alone, 

£  2 

and  minimization  of  2*  $  gives  deeper  nulls  than  does  the  minimization  of 
N  ,  n*‘ 

We  next  examine  the  case  of  two  nulls  Imposed  symmetrically  on  either  side  of 
a  null  of  the  original  pattern.  Locations  of  8  ■  15.  95*  and  17.  36*  were  chosen 
symmetrically  placed  on  either  side  of  the  null  between  the  fifth  and  sixth  sidelobes 
of  the  original  pattern  and  spaced  apart  the  same  distance  (measured  by  sin  8 )  as 
the  two  nulls  in  the  preceding  case.  Figures  10a- 10c  and  11a- He  show  the  patterns 
corresponding  to  the  perturbations  of  both  the  amplitudes  and  phases  of  *he  element 
excitations  that  minimize  respectively  the  sum  of  the  squares  of  the  relative  weight 
and  the  total  weight  perturbations.  Since  the  points  at  which  the  nulls  are  imposed 
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belong  to  adjacent  and  almost  equal  sidelobes,  and  are  approximately  symmetrical 
with  respect  to  the  null  between  the  sidelobes,  the  beam  coefficients  are  approxi¬ 
mately  equal  in  magnitude  but  opposite  in  sign  (see  Table  1).  Unlike  the  previous 
case  where  the  mainlobes  of  the  two  component  beams  of  the  cancellation  pattern 
added  in  phase  to  form  a  single  mainlobe,  here  the  cancellation  pattern  in  the  vicinity 
of  the  null  positions  is  split  into  two  lobes  of  opposite  sign.  The  fact  that  the  beam 
coefficients  are  of  opposite  sign  means  that  each  of  the  beams  must  ao  a  consider¬ 
able  amount  of  work,  as  it  were,  to  cancel  the  effect  of  the  other  beam,  instead  of 
being  helped  by  the  other  beam  as  in  the  previous  case.  As  a  result,  unlike  the  pre¬ 
vious  case,  the  beam  coefficients  are  considerably  larger  in  magnitude  than  the 
values  of  the  original  pattern  at  the  null  locations  (see  Table  1).  This  is  especially 
the  case  for  the  beam  coefficients  corresponding  to  minimization  of  the  relative 
weight  perturbations  which  are  almost  nine  times  the  magnitude  of  the  original 
pattern  values  and  three  times  as  large  as  the  beam  coefficients  corresponding  to 
minimization  of  the  total  weight  perturbations.  The  reason  why  the  minimum  rela¬ 
tive  perturbation  beam  coefficients  are  so  much  larger  than  the  minimum  total 
weight  perturbation  beam  coefficients  is  that  the  taper  of  the  mainlobe  of  the  can¬ 
celling  beams  is  much  more  gradual  in  the  former  case  than  in  the  latter.  Hence 
each  of  the  minimum  relative  weight  perturbation  beams  must  do  correspondingly 
more  work  to  cancel  the  mainlobe  component  of  the  other  beam  than  is  the  case  for 
the  minimum  total  weight  perturbation  beams. 

Note  that  even  though  the  original  pattern  at  the  null  locations  is  approximately 

N  2 

-43  dB,  the  cancellation  pattern  corresponding  to  minimization  of  Y]  (|Aw  | /a  ) 

n- 1  n  a 

rises  to  approximately  -35  dB  on  either  side  of  the  pair  of  null  locations,  and  the 
perturbed  pattern  rises  to  -33  dB.  The  location  of  the  peaks  of  the  cancellation 
pattern  away  from  the  pair  of  null  locations  (and  the  positions  of  the  peaks  of  the 
component  beams)  may  be  explained  as  follows.  Starting  at,  say,  the  left  null  posi¬ 
tion  8  *  15.  95*  and  moving  to  the  left,  the  magnitude  of  both  component  beams  de¬ 
creases,  but  the  magnitude  of  the  beam  pointed  at  0  ■  15.  95*  decreases  much  less 
rapidly  at  first  than  does  the  magnitude  of  the  beam  pointed  at  9  *  17. 36*.  Hence, 
the  two  beams  being  of  opposite  sign,  the  net  effect  is  an  increase  in  the  cancellation 
pattern.  A  similar  but  less  marked  behavior  is  also  seen  in  the  case  of  minimiza¬ 
tion  of  the  total  weight  perturbations  (compare  Figures  11a  and  11c  with  Figures 
10a  and  10c  respectively).  The  perturbation  of  the  original  pattern  in  the  vicinity 
of  the  null  positions  is  considerably  more  pronounced  for  minimization  of  the  rela¬ 
tive  weight  perturbations  than  it  is  for  minimization  of  the  total  weight  perturbations 
because  the  component  beams  of  the  cancellation  pattern  are  much  broader  in  the 
former  case  than  in  the  latter  (compare  Figure  2b  with  Figure  3b). 
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Figure  8a.  Perturbed  Pattern  With 
Two  Nulls  Imposed  at  14.54*  and 
15.  94*  With  Phase-Only  Perturbations, 

Minimizing  ltfn2.  0  ■■  -  90*  to  +90* 
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Figure  8c  .  Original  Pattern  and 
Perturbed  Pattern  With  Two  Nulls 
Imposed  at  14.54*  and  15.94*  With 
Ph-se-Only  Perturbations,  Minimizing 

vo  8.  9  ■  0*  to  25* 
n 
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Figure  10c.  Original  Pattern  and 
Perturbed  Pattern  With  Two  Nulla 
Imposed  at  15.95*  and  17.36*  With 
Combined  Amplitude  and  Phase 
Perturbations.  Minimizing  Relative 
Weight  Perturbations.  9  *  0*  to  25* 


a 


j 


40 


POWER  (OB) 


Figure  1  la.  Perturbed  Pattern  With 
Two  Nulls  Imposed  at  IS.  95  s  and 
17.  36°  With  Combined  Amplitude  and 
Phase  Perturbations,  Minimizing 
Total  Weight  Perturbations. 

9  =  -90°  to  +90° 


Figure  llh.  Cancellation  Pattern  to 
Impose  Two  Nulls  at  15.95°  and  17.36° 
Combined  Amplitude  and  Phase 
.rbations.  Minimizing  Total  Weight 
r  arbations.  9  =  -90°  to  +90° 


Figure  11c.  Original  Pattern  and 
Perturbed  Pattern  With  Two  Nulls 
Imposed  at  15.95°  and  17.36°  With 
Combined  Amplitude  and  Phase 
Perturbations,  Minimizing  Total 
Weight  Perturbations.  9  *  0°  to  25° 
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An  additional  consequence  of  the  fact  that  the  beam  coefficients  are  much  larger 
in  magnitude  than  in  the  previous  case  of  nulls  on  either  side  of  a  side  lobe  maxi¬ 
mum  is  that  the  cancellation  patterns  display  considerably  higher  sidelobes  than  in 
the  previous  rise  (compare  Figures  10b  and  lib  with  Figures  6b  and  7b  respectively). 
The  more  gradual  decrease  of  the  sidelobes  of  the  cancellation  pattern  is  especially 
noticeable  for  the  minimization  of  the  total  weight  perturbations  (compare  Figure  lib 
w<th  figure  7*»)  so  that  the  effect  on  the  Derturbed  pattern  is  seen  considerably 
further  away  from  the  null  positions  in  this  case  than  it  is  in  the  previous  case  ox 
two  nulls  imposed  symmetrically  with  respect  to  a  sidelobe  maximum  (compare 
Figure  11c  with  Figure  7  c). 

Figures  12a-l2c  and  13a-13c  show  patterns  obtained  for  the  same  two  null  loca¬ 
tions,  imposed  symmetrically  with  respect  to  a  null  of  the  original  pattern,  when 
only  the  phases  of  the  element  weights  are  allowed  to  be  perturbed.  Figures  12a- 12c 

correspond  to  minimizing  the  sum  of  the  squares  of  the  phase  perturbations,  and 

N  , 

Figures  13a- 13c  to  the  minimization  of  £  (a  $  )  .  The  qualitative  features  of 

n=  1 

these  patterns  for  8  >  0  are  very  similar  to  those  discussed  above  for  combined 
amplitude  and  phase  perturbations.  Figures  12a- 12c  corresponding  to  Figures 
10a- 10c,  and  Figures  13a- 13c  to  Figures  11a- 11c.  For  0  <  0,  of  course,  the 
phase-only  perturbation  patterns  are  dominated  by -the  pair  ot  beams  which  are 
directed  toward  the  symmetric  positions  of  the  null  locations,  and  which  have  signs 
opposite  to  the  corresponding  cancellation  beams  directed  at  the  null  locations 
themselves. 

Finally  we  examine  the  case  of  two  nulls  imposed  asymmetrically  within  a  side- 
lobe.  Locations  of  9  »  IS.  23*  and  IS.  94*  were  chosen  corresponding  respectively 
to  the  peak  and  right  3  dB  point  of  the  fifth  sidelobe  of  the  unperturbed  pattern. 
Figures  14a- 14c  and  15a- 15c  show  the  patterns  for  combined  amplitude  and  phase 
perturbations  that  minimize  respectively  the  relative  and  total  weight  perturbations. 

In  contrast  to  the  case  of  two  nulls  imposed  symmetrically  with  respect  to  a  sidelobe 
maximum,  here  the  mainlobe  of  the  cancellation  pattern  for  minimized  relative 
weight  perturbations  is  seen  to  be  split  into  two  lobes  (compare  Figure  14b  with 
Figure  6b).  Referring  to  Table  1  we  see  that  the  coefficients  of  the  two  beams  com¬ 
posing  the  cancellation  pattern  are  of  opposite  sign  and  have  magnitudes  more  than 
five  times  larger  than  the  corresponding  unperturbed  pattern  values,  indicating  that 
the  two  beams  are  doing  considerable  work  in  cancelling  the  effect  of  each  other  in 
addition  to  cancelling  the  original  pattern.  This  behavior  of  the  cancellation  pattern 
may  be  explained  qualitatively  as  follows.  The  mainlobe  of  the  beam  directed  at 
the  location  of  the  sidelobe  peak  at  IS.  23*  is  much  broader  than  the  sidelobe  of  the 
original  pattern  and  hence  tapers  off  much  more  gradually  than  does  the  unperturbed 
sidelobe.  Hence,  if  the  beam  directed  at  15.  23*  was  exactly  matched  in  magnitude 
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to  the  original  pattern  there  it  would  have  a  component  considerably  larger  in  mag¬ 
nitude  than  the  value  of  the  unperturbed  pattern  at  IS.  94*.  To  counter  this  unwanted 
"overshoot",  as  it  were,  of  the  beam  directed  at  15.  23*,  the  beam  directed  at 
IS.  94*  must  be  of  opposite  sign.  The  component  at  15. 23*  of  the  mainlobe  of  the 
beam  directed  at  15.  94*  in  turn  requires  cancellation  by  the  peak  of  the  beam  dir¬ 
ected  at  15. 23*  which  forces  the  magnitude  of  the  latter  beam  to  exceed  the  magni¬ 
tude  of  the  original  pattern  at  15.  23*. 

A  similar  but  less  pronounced  behavior  of  the  cancellation  pattern  occurs  for 
minimized  total  weight  perturbations  (compare  Figure  15b  with  Figure  7b).  The 
mainlobe  of  the  beam  directed  at  15.  23*  in  this  case,  being  less  wide  than  the  main- 
lobe  of  the  minimized  relative  weight  perturbation  beam  but  still  twice  as  wide  as 
the  sidelobe  of  the  original  pattern,  has  a  smaller  "overshoot"  at  15. 9**  requiring 
correspondingly  less  compensation  by  the  beam  directed  at  15.  94*. 

Figures  16a-  16c  and  17a- 17c  show  the  patterns  obtained  for  the  case  of  the  same 


two  null  locations  at  15.  23*  and  15.  94*  when  only  the  phases  of  the  weights  are 

N  „ 

allowed  to  vary.  Figures  16a-16c  correspond  to  minimization  of  £  p  and 

N  2  n=l  n 

Figures  17a- 17  c  to  the  minimization  of  V  (a  p  ).  The  behavior  of  these  patterns 

n*  1  n  n 

for  0  >  0  is  similar  to  that  of  the  corresponding  combined  amplitude  and  phase 


perturbation  patterns. 


Following  this  discussion  of  the  basic  characteristics  of  nulling  at  selected 
points,  we  now  discuss  a  set  of  computations  performed  to  investigate  cancellation 
in  a  given  sector  of  a  pattern.  We  define  the  power  cancellation  ratio  in  the  inter- 
valAMfjiJS  «a 


where  pQ  <0 )  is  the  original  pattern  and  p<0 )  is  the  perturbed  pattern. 
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Figure  14a.  Perturbed  Pattern  With 
Two  Nulls  Imposed  at  15.23*  and 
15.  94*  With  Combined  Amplitude  and 
Phase  Perturbations,  Minimizing 
Relative  Weight  Perturbations. 

8  =  -90*  to  +90* 


Figure  14b.  Cancellation  Pattern  to  B 

Impose  Two  Nulls  at  15.  23*  and  IS.  94*  8 

With  Combined  Amplitude  and  Phase  m 

Perturbations,  Minimizing  Relative  Weight  5 
Perturbations.  8  »  -90*  to +90*  ® 


Figure  14c.  Original  Pattern  and 
Perturbed  Pattern  With  Two  Nulls 
Imposed  at  15.  23*  and  15.94*  With 
Combined  Amplitude  and  Phase 
Perturbations,  Minimizing  Relative 
Weight  Perturbations.  8  *  0*  to  25* 
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Figure  15c.  Original  Pattern  and 
Perturbed  Pattern  With  Two  Nulls 
Imposed  at  15.23*  and  15.94*  With 
Combined  Amplitude  and  Phase 
Perturbations,  Minimizing  Total 
Weight  Perturbations.  9  *  0*  to  25* 
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The  procedure  used  to  study  the  behavior  of  the  cancellation  ratio  for  the  four 
types  of  weight  perturbations  was  as  follows.  Starting  with  the  41-element,  40  dB 
Chebyshev  taper  pattern,  the  left  hand  boundary  of  the  sector  of  interest,  8  j,  was 
fixed  at  the  peak  of  the  fifth  sidelobe  at  6  «  IS.  23*.  A  second  location  6  j.  was 
then  found  such  that  the  cancellation  ratio  in  the  interval  [15.  23*,  6  jl  was  30  dB 

when  a  null  was  imposed  at  IS.  23*  and  8  ,  using  combined  amplitude  and  phase 

N  ,2 

perturbations  and  requiring  J'  |AwJ  to  be  minimized.  This  gave  fl,  -  15.78*. 

n«  1  n  2 

The  cancellation  ratio  for  all  four  types  of  weight  perturbations  in  the  sector 

(IS.  23*,  15.78*]  with  nulla  placed  at  15. 23*  and  15. 78*  was  then  determined.  Next 
the  sector  was  enlarged  to  the  right  by  adding  a  length  (measured  by  sin  0 )  equal 
to  sin  (15.78*)  -  sin  (15. 23*)  giving  8  2  *  16.  33*.  The  cancellation  ratio  in  the 
sector  [15.  23*,  16. 33*]  was  then  calculated  for  all  four  types  of  perturbations  with 
nulls  constrained  to  lie  at  15.  23*.  15.  7S*.  and  16.33*.  The  procedure  was  then 
continued  .'u  the  same  way  by  adding  equispaced  nulls.  The  cancellation  ratios  ob¬ 
tained  are  shown  in  Table  2.  Table  3  shows  the  values  of  beam  coefficients  and  null 
depths  for  these  cases. 

Referring  to  Table  2,  we  see  that  the  cancellation  ratios  obtained  with  the  two 
types  of  combined  amplitude  and  phase  perturbations  decrease  steadily  as  the  number 

of  nulls  increases  indicating  more  effective  cancellation.  Slightly  lower  ratios  are 

N  N 

obtained  by  minimizing  |  A w  I  2  than  by  minimizing  Y\  (|  Aw  | /a  )2.  The 

n*l  n  n  =  i  n  n 

improvement  (that  is,  decrease)  of  the  cancellation  ratio  with  the  increase  in  number 

of  nulls  may  be  explained  by  noting  that  in  requiring  nulls  to  be  placed  at  various 
locations  we  are,  equivalently,  requiring  the  cancellation  pattern  to  exactly  match 
the  original  pattern  at  these  locations.  As  the  number  of  points  at  which  the  cancella¬ 
tion  pattern  is  matched  to  the  original  pattern  increases,  the  better  is  the  overall 
fit  of  the  cancellation  pattern  to  the  original  pattern  and  hence  the  cancellation  ratio 
decreases,  indeed,  in  the  limit  of  using  all  the  degrees  of  freedom  et  our  disposal 
by  imposing  nulls  at  N  locations,  the  cancellation  pattern  will  exactly  match  the 

original  pattern  everywhere  and  hence  the  cancellation  ratio  will  be  minus  infinity. 

N  , 

The  slightly  lower  Cancellation  ratios  obtained  by  minimizing  T'  I  Aw  I c  as  com- 

N  ,  ns  1  ™ 

pared  with  those  obtained  by  minimizing  (|  Aw  I  /a)  is  consistent  with  what  we 
1  n=  1  n  n 

have  observed  above’  in  our  examination  of  the  patterns  corresponding  to  nulling  at 
one  or  two  selected  points.  The  cancellation  beams  used  when  the  relative  weight 
perturbations  are  minimized  have  a  much  broader  mainlobe  than  those  used  when 
the  total  weight  perturbations  are  minimized,  and  hence  are  more  likely  to  result 
in  higher  perturbed  pattern  sidelobe  values  in  between  the  null  locations. 
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Table  3.  Values  of  the  Imposed  Null  Locations  <0nuu>.  Original  Pattern  (pQ).  Beam  Coefficients  (b).  and 
Perturbed  Pattern  (p)  for  Four  Types  of  Weight  Perturbations.  (Numbers  in  parentheses  are  powers  of  ten) 


The  behavior  of  the  cancellation  ratios  obtained  for  the  two  types  of  phase-only 

perturbations  is  in  marked  contrast  to  the  behavior  of  the  cancellation  ratios  for 

the  two  types  of  combined  amplitude  and  phase  perturbations.  The  cancellation 

N  , 

ratios  obtained  by  minimizing  T\  0  increase  as  the  number  of  null  locations  in¬ 
ns  1  n  N  2 

creases  beyond  three,  and  those  obtained  by  minimizing  T'  (a  p  )  increase  as 

n=  1  n  n 

the  number  of  null  locations  increases  beyond  two.  For  six  null  locations  there  is 
virtually  no  effective  cancellation  by  either  scheme  of  phase-only  perturbations. 

In  general  the  cancellation  ratios  obtained  by  minimizing  the  sum  of  the  squares  of 
the  phase  perturbations  themselves  are  lower  than  those  obtained  by  minimizing 
the  sum  of  the  squares  of  the  weighted  phase  perturbations.  The  reason  for  the 
much  poorer  cancellation  obtained  with  the  phase-only  perturb.' •  x  schemes  as 
compared  with  the  combined  amplitude  and  phase  perturbation  schemes  is  the  small 
angle  approximation  used  to  linearize  the  phase-only  nulling  problem.  As  the  number 
of  imposed  nulls  increases  and  coupling  between  the  component  beams  of  the  phase- 
only  cancellation  patterns  increases,  the  beam  coefficients  become  larger  and  larger 
in  magnitude  (see  Table  3)  and  consequently  the  phase  perturbations  increase  in 

magnitude.  For  example,  for  two  nulls  at  15.23°  and  15.78°,  the  largest  phase 

N  2  N  2 

perturbations  obtained  by  minimizir.'  £  0n  and  V  <an0n'  are  1.  55°  and  5.47° 

respectively.  In  contrast,  for  six  nulls  the  largest  phase  perturbation  is  42. 7  °  and 

N  2 

18  phase  perturbations  are  greater  than  10*  for  minimizing  T'  p  while  for  mini- 

N  2  n=  1  n 

mixing  (a  d  )  the  largest  phase  perturbation  is  65°  and  10  phase  perturba- 
n=  1 

tions  are  greater  than  10°.  The  larger  the  phase  perturbations,  the  worse  the  small 
angle  approximation.  Hence  the  null  depth  achieved  degrades  (see  Table  3)  and 
consequently  the  overall  sector  cancellation  is  less  effective.  The  fact  that  the 
small  angle  approximation  becomes  less  valid  as  the  phase  perturbations  grow,  also 

explains  why  somewhat  better  cancellation  ratios  are  obtained  by  minimizing 
N  2  N  2 

£  p*  than  by  minimizing  £  (an  d  ).  The  phase  perturbations  obtained  with  the 
n* 1  n«  1 

former  scheme  are,  in  a  mean  squares  sense,  lower  than  those  obtained  with  the 
latter  scheme  and  so  satisfy  the  small  angle  approximation  better. 


4.  GENERALIZATIONS 


In  Section  2  we  saw  that  the  process  of  determining  weight  perturbations  to 
place  nulls  at  a  set  of  specified  locations  led  to  a  representation  of  the  cancellation 
pattern  in  terms  of  beams.  Minimization  of  the  sum  of  the  squares  of  the  weight 
perturbations  relative  to  the  original  weights  led  to  a  representation  of  the 
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cancellation  pattern  in  terms  of  beams  corresponding  to  an  element  amplitude  taper 
equal  to  the  square  of  the  original  taper,  while  minimization  of  the  sum  of  the 
squares  of  the  total  weight  perturbations  led  to  a  representation  in  terms  of  beams 
corresponding  to  a  uniform  amplitude  taper.  We  also  saw  that  the  latter  minimiza¬ 
tion  problem  is  equivalent  to  minimizing  the  mean  square  difference  of  the  can¬ 
cellation  pattern  over  a  full  period  of  the  variable  u  »  kd  sin  (0 ).  It  is  natural  to 
inquire  whether  there  is  a  minimization  of  some  property  of  the  cancellation  pattern 
that  corresponds  to  the  minimization  of  the  sum  of  the  squares  of  the  relative  weight 
perturbations.  The  object  of  this  section  is  to  answer  this  question  and  in  so  doing 
to  generalize  some  of  the  results  obtained  in  Section  2. 

We  begin  by  returning  to  the  problem  addressed  in  Section  2  of  determining 
combined  amplitude  and  phase  perturbations  to  place  nulls  at  specified  locations  in 
a  given  pattern.  Noting  the  general  relation  between  the  weights,  wn>  of  the 
elements  of  an  equispaced  linear  array  and  the  pattern,  p(u)  of  the  array 


N  f 


where 


N 


w  e  , 

n=  1  n 

dn  =  ~2~  —  U 


and 

u  =  kdsin  6 

it  follows  that  the  sum  of  the  squares  of  the  total  perturbations  Awq  equals  the  mean 

of  the  squared  cancellation  pattern,  hp(u)  «  p(u)  -  pQ(u);  that  is,  _ _ 

N  n 

IhwJ2  =  J'  |p(u)  -  pQ(u)|2  du  .  (28) 

n_  ~n 

Hence  minimizing  the  sum  of  the  squares  oi  the  total  weight  perturbations  is  equiva¬ 
lent  to  minimizing  the  mean  squared  cancellation  pattern. 

We  now  ask  whether  there  is  a  similar  pattern  minimization  that  corresponds 
to  the  minimization  of  the  sum  of  the  squares  of  the  weight  perturbations  relative 
to  the  original  weights: 
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2  N 


£.-£*-  -  £,(£*!)  ■ 
a=  1  _  ^  n  s  n=  1  '  n  ' 


n=  1  .  J  n  s 

a  e 
n 


In  fact  the-e  is,  and  the  answer  is  given  in  terms  of  minimizing  the  mean  square 
of  the  convolution  of  the  cancellation  pattern  with  another  pattern.  For  in  general, 
if 

N  ■  , 

-id„u 


V1  “Jun 

Pj<u)  =  2~!  wln  e 


A  -jd  u 

P2(U)  =  2^  w2n  e  n 

n=  l 


then  it  is  easy  to  verify  by  direct  substitution  that  the  convolution  of  p^  and  p ^ 
defined  by 

7r  rr 

tPl  *P2]  (u)  = Pj(u  )  p2(u  -  u  )  du  =  2^  J*  Pj  (u  -  u’) p2 <u’) du’ 


has  the  representation 

N  -id  u 

[Pl*P2l(u)=  E  wlnw2ne  *  n  . 


Hence,  from  Eq.  (27), 

7T  . 


2?  J  |lPl*P2)  (u)|  d“  =  ^  |wlnW2n>: 


It  follows  by  letting  w.  =  Aw  and  wn  =  1/a  that 
in  n  2n  n 

n?i(^T"^=^  f  ll<p-po)*gl(u>|  du 


where 


g(u)  =  E  T- 


'  a 

n=  1  n 


1  /^n11 

—  e 
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Thus  the  minimization  of  the  sum  of  the  squares  of  the  relative  weight  perturbation 
is  equivalent  to  the  minimization  of  the  mean  square  of  the  convolution  of  the  can¬ 
cellation  pattern  with  the  pattern  corresponding  to  element  weights  equal  to  the 
reciprocals  of  the  amplitudes  of  the  original  weights. 

Note  that  we  can  also  recast  the  equivalence  of  minimizing  the  3um  of  the 
squares  of  the  total  weight  perturbations  with  minimizing  the  mean  squared  can¬ 
cellation  pattern  [Eq.  (38)]  by  letting  wJn  -  iw^  and  =  1  in  Eq.  (29)  thereby 
obtaining 


II 

H  lAwnl2^2^  f  i  HP  -  PQ> g]  (u) I  2  du 
n=  1  ‘  J 


V  ‘j<1nU 
g(u)  =  2-j  e 
n=  1 


I  Nu\ 
\  2  ) 

sin  (  j 


Hence,  minimizing  the  sum  of  the  squares  of  the  total  weight  perturbations  is 
equivalent  to  minimizing  the  mean  square  of  the  convolution  of  the  cancellation 
pattern  with  sin  (  j  /sin  (  y  )  •  The  result,  Eq.  (30),  can  also  be  obtained  di¬ 
rectly  from  Eq.  (28)  by  noting  (as  is  easily  verified  using  Eq.  (31)]  that  the  convolution 
of  an  cquispaced  linear  array  pattern  with  the  function  sin  /sin(  ^  )  is  the 

array  pattern  itself  so  that 
t  Nu  \ 

SUM  ~T  j 

(p  -  p  )  * -  (u)  =  p(u)  -  p  (u)  . 

L  ° 


(u)  =  p(u)  -  PQ(u)  . 


We  can  now  go  further  and  observe  that  the  two  sums  of  squares  of  weight  perturba¬ 
tions, 

El-/  -  E,  (A1) 

whose  minimization  in  conjunction  with  placing  nulls  at  a  set  of  specified  locations 
led  to  the  two  forms  of  cancelling  beams  treated  in  this  report,  are  special  cases 
of  a  more  general  weighted  sum  of  squares 


J-i  /|  Aw  |\“  1  r  1  I  2 

2-.  1 — - — )  =  77  j  l|(p  -  po>  :;g]  (u)r  du 

n=  1  \  n  /  J 
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where 


«<«)«  £  T-  e‘JdnU 
n  =  1  n 


ft  is  then  reasonable  to  ask  whether,  given  any  set  of  non-zero  real  coefficients 

c  we  can  find  the  set  of  perturbations  Aw  that  will  result  in  the  nulls  at  a  set  of 

n  ..  v  n 

prescribed  locations  u » u^,  k*=l,2,  ....  M,  and  that  will  minimize 
N 


£  (W 


or,  equivalently,  minimize 

»  2 
J*  ![<P-P0)*g](u)l  du 


•it 


The  condition  for  nulls  in  the  perturbed  pattern  at  u*  u^,  k=  1,  2,  . . , ,  M,  is 
0  *  pit^)  =  P0(uk)  +  Ap(uk) 

V1  "J  dn  uk 
“  po<V  +  nt-,  Awn  6 


or 


N 


2d  e  K  Aw  *  -  PjjiUjj).  k  *  U  2 . M 

n*  1 


(33) 


so  that  we  wish  to  find  the  set  of  weight  perturbations  Awn  satisfying  Eqs.  (33)  and 
minimizing 

JC  (Ia wnl /cn)2  . 

The  problem  can  be  expre  „ed  in  the  form  of  finding  the  vector  of  weight  perturba¬ 
tions 


Aw  *  [AWj,  Aw2,  ....  iwN| 
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satisfying 


A  Aw  ■  jr 

and  minimising 


T 


(Aw)1  C  Aw 


where  A  is  the  HXN  matrix 
*lul 


Jd.u,  -Jd2u 

e  ...  e 


/JdauM 

6  C 


C  is  the  N  X  N  diagonal  matrix 
.  -2 


r2  0 


0 


-2 


N 


and 


-JVi 
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(34) 


10, 


Using  a  result  from  the  theory  of  generalized  inverses,  the  solution  is  given  by 
Aw  «  C*1  A*  (AC*1  aV1  i 

so  that  Aw  is  a  'inear  combination  of  the  columns  of  C* 1 A^  with  coefficients  from 
the  vector 


b  »  (AC^A*)*1* 


But 


C_iA^ 


cjJVl 


c2^2*1 


cNe 


JVi 


idi  ua 


2  J“l “2 
Cj  e 


2  Jd2u2 
c2  e 
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2  J“1 “M 
Cj  e 


c2eld2“" 


2  cidNUM 
CN  C 


so  that 


M 


Aw_  *  c2  J]  b__  e 
m=  1 


jdu 
n  m 


and  the  perturbed  pattern 


p(u)  -  po(u)  +  A  p(u) 


where 


M 


N 


APCU,  ■  Zb  7 


_  C  e 

m-  1  m  n«  1  n 


2  -Jd„<«-Um) 


If  the  original  pattern  pQ(u)  is  real  as  we  assumed  in  Section  2,  then  the  cancellation 
pattern  must  likewise  be  real.  This  is  assured  by  assuming  even  symmetry  of  the 
weighting  coefficients  cR  with  respect  to  the  phase  reference  center;  that  is. 


cn  "  ®N  -  n+  l*  n  *  U  2* - *  N 


Then  the  cancellation  pattern  is  represented  as  the  sum  of  M  real  beams,  one  pointed 

at  each  of  the  null  locations,  corresponding  to  an  amplitude  taper  equal  to  the  square 

of  the  weighting  coefficients  cQ.  The  two  types  of  beams  treated  in  Section  2  are 

special  cases  of  this  result.  Minimum  relative  weight  perturbations  are  obtained, 

2  2 

as  noted  above,  by  letting  cn  «  aR  which  gives  cn  ■  aR  and  a  representation  of 
the  cancellation  pattern  in  terms  of  beams  corresponding  to  an  amplitude  taper  equal 
to  the  square  of  the  original  taper.  Minimum  total  weight  perturbations  are  obtained 
by  letting  cR  »  1  so  that  cR  ■  1  and  the  cancellation  pattern  is  represented  as  the 
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sum  of  beams  corresponding  to  a  uniform  amplitude  distribution  (that  is, 

sin  (  ^  )  /sin(  i  )  beams).  We  could  equally  well  consider  other  types  of  beams. 

*  “  aw  2 

For  example,  by  letting  c  =  va  ,  we  obtain  c  *  a_  so  that  the  cancellation 
beams  are  replicas  of  the  original  pattern  itself. 

/  decision  regarding  which  type  of  beams  to  use  for  nulling  in  a  particular 
situation  must  take  into  account  the  characteristics  of  the  cancelling  beams  and 
how  they  affect  the  perturbations  of  the  original  pattern  at  points  other  than  the 
imposed  null  locations.  We  saw  in  Section  3  that  the  deviations  of  the  perturbed 
pattern  from  the  original  pattern  at  points  other  than  the  null  locations  depend 
strongly  on  the  width  of  the  mainlobe  and  the  taper  of  the  sidelobes  of  the  cancelling 
beams.  Because  of  the  inverse  relationship  between  width  of  mainlobe  and  taper  of 
sidelobes,  there  is  a  tradeoff  between  relatively  small  perturbations  spread  out  for 
a  relatively  large  distance  from  the  null  locations  (as  with  the  sin  )  /sin  (-fa) 
beams),  and  relatively  large  perturbations  restricted  to  a  relatively  small  vicinity 
of  the  null  locations  (as  with  the  beams  corresponding  to  a  taper  proportional  to  the 
square  of  an  already  highly  tapered  element  amplitude  distribution).  It  is  possible 
that  the  relationship  established  in  this  section  between  the  type  of  cancelling  beams 
and  the  minimization  of  the  convolution  of  the  cancellation  pattern  with  a  pattern 
uniquely  associated  with  the  type  of  cancelling  beams,  can  be  used  to  help  clarify 
and  quantify  this  tradeoff  between  weak,  distributed  perturbations  and  strong,  localized 
perturbations,  but  more  work  remains  to  be  done  here. 

Although  we  have  so  far  in  this  section  restricted  our  attention  to  combined 
amplitude  and  phase  perturbations,  the  results  established  here  can  be  extended  in 
a  parallel  way  to  phase-only  nulling  as  well.  Equation  (32)  can  be  used  as  the  start¬ 
ing  point  for  an  analysis  of  phase-only  nulling,  with  the  form  of  the  weight  perturba¬ 
tions  given  by  | 


Aw  aw  —  w 
n  n  on 


w  e  -  w 
on  on 


(e  "  -  l) 


where  0n  is  the  phase  perturbation  of  the  nxn  weight.  For  small  perturbations. 


Aw  *  'I*  6 

n  4  onyn 
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special  cases  treated  in  Section  2,  the  minimization  of  and 

correspond  to  the  choice  of  cR  *  an  and  cn  *  1  respectively,  for  which  Eq.  (37) 
yields  Eq.  (19)  and  Eq.  (23). 

The  cancellation  pattern  corresponding  to  Eq.  (37)  is 

£  "^dn u 
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}• 


For  if/  -  d  u  ,  this  gives  the  approximate  representation  of  the  phase-only  cancella¬ 
tion  pattern  as  the  superposition  of  M  pairs  of  beams,  one  member  of  the  pair 
directed  at  a  null  location,  and  the  other  member,  of  opposite  sign,  directed  at  a 
location  symmetric  with  respect  to  the  mainlobe  of  the  original  pattern. 


3.  CONCLUSIONS 

In  this  report,  we  have  considered  the  problem  of  imposing  nulls  in  the  pattern 
of  a  linear  array  of  equispaced,  isotropic  radiators  subject  to  the  condition  that  the 
perturbations  of  the  element  weights  be  minimized.  In  Section  2  we  analyzed  four 
forma  of  the  problem:  (la)  nulls  imposed  with  combined  amplitude  and  phase  per¬ 
turbations,  minimizing  the  sum  of  the  squares  of  the  perturbations  relative  to  the 
original  weights;  (lb)  nulls  imposed  with  combined  amplitude  and  phase  perturbations, 
minimizing  the  sum  of  the  squares  of  the  total  weight  perturbations;  (2a)  nulls  im¬ 
posed  with  phase  perturbati.ns  only,  minimizing  the  sum  of  the  squares  of  the  per¬ 
turbations  relative  to  the  original  weights;  and  (2b)  nulls  imposed  with  phase 
perturbations  only,  minimizing  the  sum  of  the  squares  of  the  total  weight  pertuiba- 
tions.  The  results  are  summarized  in  Table  4.  Forms  (la)  and  (lb)  can  be  solved 
exactly  for  the  required  weight  perturbations.  Form  (la)  leads  to  a  representation 
of  the  cancellation  pattern  as  a  superposition  of  beams,  one  directed  at  each  im¬ 
posed  null  location,  each  beam  corresponding  to  an  element  taper  proportional  to 
the  square  of  the  original  kaper;  while  form  (lb)  leads  to  a  representation  of  the 
cancellation  pattern  as  a  superposition  of  beams  of  the  sin  (  )  /sin(  •j  )  typo. 
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Table  4.  Summa 


Forms  (2a)  and  (2b)  cannot  be  solved  exactly,  but  approximate  solutions  can  be 
obtained  by  assuming  that  the  required  phase  perturbations  are  small.  Minimiza¬ 
tion  of  the  relative  weight  perturbations  is  then  approximated  by  the  minimization 
of  the  sum  of  the  squares  of  the  phase  perturbations,  while  minimizing  the  total 
weight  perturbations  is  approximated  by  minimizing  the  sum  of  the  squares  of  the 
products  of  the  phase  perturbations  with  the  amplitudes  of  the  weights.  The  result¬ 
ing  cancellation  patterns  can  then  be  approximately  represented  as  the  sum  of  pairs 
of  beams,  one  pair  for  each  null  location.  One  member  of  the  pair  of  beams  is 
directed  at  the  null  location,  while  the  other,  of  opposite  sign,  is  directed  at  the 
location  symmetric  with  respect  to  the  direction  of  the  mainlobe  axis  of  the  original 
pattern.  The  form  of  beams  is  the  same  as  that  for  the  respective  combined  ampli¬ 
tude  and  phase  perturbation  problem;  that  is,  beams  corresponding  to  a  taper  equal 
to  the  square  of  the  original  amplitude  taper  for  (2a)  and  sin  (  )  /sin^-j  j  beams 

for  (2b). 

In  Section  3  we  presented  and  discussed  numerical  calculations  performed  to 
display  the  principal  features  of  the  solutions  obtained  in  Section  2.  The  original 
pattern  was  taken  to  be  that  of  a  41  element  array  with  half  wavelength  spacing  and 
a  40  dB  Chebyshev  taper.  We  examined  tiie  patterns  obtained  with  the  four  methods 
of  Section  2  for  the  cases  of  one  imposed  null  at  the  peak  of  a  sidelobe,  two  nulls 
imposed  symmetrically  at  the  -3  dB  points  of  a  sidelobe,  two  nulls  imposed 
symmetrically  on  either  side  of  a  null  of  the  original  pattern,  and  two  nulls  imposed 
asymmetrically,  one  at  the  peak  of  a  sidelobe  and  the  other  at  a  -3  dB  point  of  the 
sidelobe.  We  saw  that  in  general  the  "tapered"  beams  corresponding  to  forms  (la) 
and  (2a)  resulted  in  larger  perturbations  of  the  original  pattern  in  the  vicinity  of  the 
Imposed  null  locations  and  smaller  perturbations  far  away  from  the  imposed  nulls 
than  did  the  "sine"  beams  corresponding  to  forms  (lb)  and  (2b).  This  was  a  conse¬ 
quence  of  the  fact  that  the  "tapered"  beams  had  wider  mainlobes  and  lower  side- 
lobes  than  did  the  "sine"  beams.  Nulling  in  the  phase-only  patterns  was  accom¬ 
panied  by  an  approximately  equal  and  opposite  distortion  of  the  original  pattern  at 
the  points  symmetric  to  the  imposed  null  locations  with  respect  to  the  axis  of  the 
original  pattern  mainlobe.  As  a  result,  both  amplitude'and  phase  weighting  are 
required  to  produce  nulls  that  are  symmetric  with  respect  to  the  direction  of  the 
mainlobe.  - 

In  Section  3  we  also  discussed  the  results  of  calculations  performed  to  examine 
cancellation  within  a  sector  of  the  pattern  as  the  width  of  the  sector  was  increased 
by  adding  equispaced  imposed  nulls.  We  saw  that  the  cancellation  achieved  by  the 
two  combined-amplitude-and-phase  perturbation  methods  became  increasingly 
effective  as  the  width  of  the  sector  was  increased,  but  that  the  depth  of  nulls  and 
cancellation  achieved  by  the  two  phase-only  methods  deteriorated  severely  as  the 
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width  of  the  sector  was  increased.  This  deterioration  of  the  performance  of  the 
phase-only  methods  was  attributed  to  the  fact  that  as  imposed  nulls  were  added, 
interference  between  the  cancelling  beams  became  more  pronounced,  the  required 
phase  perturbations  were  larger,  and  consequently  the  small  angle  approximation 
used  to  derive  the  phase-only  solutions  became  more  and  more  inaccurate 

In  Section  4  we  considered  the  problem  of  finding  a  pattern  minimization  \'hat 
was  equivalent  to  minimizing  the  relative  weight  perturbations,  in  the  same  way 
that  minimizing  the  mean  square  cancellation  pattern  was  equivalent  to  minimizing 
the  total  weight  perturbations.  We  found  that  there  was  an  equivalent  pattern 
minimization,  namely,  minimizing  the  mean  square  of  the  convolution  of  the  can¬ 
cellation  pattern  with  the  pattern  corresponding  to  element  weights  equal  to  the 
reciprocals  of  the  original  amplitudes.  We  then  showed  that  the  types  of  minimiza¬ 
tion  considered  in  Section  2  w»re  special  cases  of  a  more  general  minimization  of 
the  sum  of  squares  of  the  weigh'  perturbations  divided  by  arbitrary  non-zero  real 
coefficients— a  minimization  whUh  had  its  equivalent  counterpart  in  pattern  space 
in  terms  of  minimizing  the  mean  square  of  the  convolution  of  the  cancellation 
pattern  with  the  pattern  corresponding  to  element  weights  equal  to  the  dividing 
coefficients.  We  then  solved  this  generalized  minimization  problem  for  combined 
amplitude  and  phase  perturbations,  and  for  phase-only  perturbations. 
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Appendix  A 

Nulling  With  Minimum  Praam  Perturbation* 


In  this  appendix  we  solve  the  following  problem: 

Given  a  linear  a '•ray  of  N  equispaced  isotropic  radiators  and  an  initial  set  of  com¬ 
plex  element  weights,  w^,  n=  1.  2 . N.  determine  the  set  of  complex  weights 

wn,  n=  1,  2,  ....  N.  which  gives  the  array  pattern  with  nulls  at  a  set  of  M  pres¬ 
cribed  pattern  locations,  u^,  k  =  1.2.....  M,  M  <  N,  and  which  differs  as  little  as 
possible  in  the  mean  sense  from  the  initial  pattern.  When  we  use  the  word  "mean" 
here  we  define  the  interval  over  which  the  mean  Is  taken  to  be  the  visible  region 
-  s  e  -  ^  so  that  the  variable  u  =  kd  sin  0  is  to  be  integrated  from  -kd  to  +kd. 

The  initial  pattern  pQ(u)  is  given  in  Section  2,  as 


(u)  =  T  w„ 


and  the  perturbed  pattern  p(u)  by 


p<u)  =  52 


The  difference  between  the  patterns  is  then 


E-J  a  i 

.  (wn  ’  won)  e 
n  —  l 


( 


and  the  mean  square  difference  between  the  patterns  is 
kd  N  N 

|p(u)  -  p0<u)  | 2  du  -  E  E  -  w*m)  <wn  -  won) 
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where 


^  "  «wl  -  wol'  W2  *  wo2*  —  WN  *  woNl 

,T 


(Aw,.  Aw, 


awn( 


and 


[C  ]  Bin  '  •lnc«dn-dm)“J-  [  C] 


tun 


Note  that  for  half  wavelength  spacing,  kd  *  s  and  C  becomes  the  identity  matrix  so 
that 

kd 

TET  S  IpM  ~  PQtu)l 2  *  (Aw)t  (Aw) 


-kd 
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-  E 
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n=  1 


(Al) 


The  perturbed  pattern  is  required  to  have  nulls  at  the  M  locations  u  =  u^, 
k  »  1,  2,  ....  M.  so  that 
N  .  . 
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or 


(*idlUk.  . e^  ukj  it  =  1,  2 . M 


so  that 


*n  *  bk 


Jdn“k 


This  expression  for  the  weight  perturbations  is  identical  to  Eq.  (10)  obtained  by 

requiring  £|  Awn|  2  to  be  a  minimum.  Of  course  this  conclusion  is  also  an 
immediate  consequence  of  Eq.  (Al).  but  it  is  of  interest  to  obtain  it  by  solving  the 
minimum  pattern  perturbation  problem  with  arbitrary  uniform  element  spacing. 

In  the  case  of  phase-only  element  control,  half  wavelength  spacing,  and  small 
phase  perturbations,  starting  with  Eq.  (Al)  we  have 


I  p(u)  -  P0<u)l 1  du 


E  I  Aw  |2 
n=  1  n 

N 


so  that  minimizing  £(an  $n’  ls  equivalent  to  minimizing  the  mean  square  pattern 
perturbation. 

j  In  the  above,  as  mentioned  earlier,  we  have  defined  the  mean  with  respect  to 
the  visible  region  of  the  array  pattern.  Since  the  array  pattern  for  an  equispaced 
linear  array  is  periodic  with  period  2  s  in  the  variable  u  =  kd  sin  6  ,  it  is  also 
possible  to  define  the  mean  with  respect  to  one  period  of  the  variable  u  regardless 
of  whether  or  not  the  visible  region  constitutes  one  complete  period.  If  this  is  done, 
then  the  mean  square  difference  between  the  original  and  perturbed  pattern  is 

i  S  lp(U)  '  Po(u)l  2  dU  *  E^  (wm  -wom)  (wn  “  w«>  f  e  ^  dm>U  1 

■f  ~  -» 

N  N 
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so  that  the  sum  of  the  squares  of  the  weight  perturbations  equals  the  me;  n  square 
of  the  cancellation  pattern,  the  mean  being  taken  with  respect  to  a  full  pe  riod  of  the 
variable  u  =  kd  sinO  .  Thus  a  least  mean  square  match  to  an  initial  element  ex¬ 
citation  is  equivalent  to  a  least  mean  square  match  over  one  period  of  the  initial 
pattern.  In  the  special  case  of  half  wavelength  spacing,  d  *  A/ 2,  the  period 
-  n  s  kd  sin  0  -  it  coincides  with  the  visible  region  and  Eq.  (A4)  becomes  identical 
with  Eq.  (Al).  When  d  >  A/2,  then  a  match  over  a  full  period  of  the  variable  u  im¬ 
plies  a  match  over  an  angular  sector  smaller  than  the  visible  region.  This  how¬ 
ever,  is  not  too  serious  since  the  pattern  is  periodic  outside  this  sector. 
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